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In the long-wavelength limit fc = 0, the response function has been investigated with respect to the 
external and internal fields which is expressed by the external and internal conductivity, respectively. 
Molecular dynamics (MD) simulations are performed to obtain the current-current correlation func- 
tion and the dynamical collision frequency which are compared with analytical expressions. Special 
attention is given to the dynamical collision frequency and the description of plasma oscillations 
in the case of fc = 0. The relation between the external and internal conductivity and to the 
current-current correlation function is analyzed. 
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I. INTRODUCTION 

The treatment of strongly correlated Coulomb systems is a challenge for many-particle theories. It has applications 
in different fields such as dense ionic plasmas and the electron-hole plasma in excited semiconductors. Within a 
quantum statistical approach, the methods of equilibrium and non-equilibrium Green functions have successfully 
been utilized to calculate the properties of dense plasmas, see Q. However, a problem is the validity of perturbative 
approximations when using the Green function approach for strongly correlated systems. 

With increasing computer capacities, simulation techniques such as molecular dynamics (MD) simulations have 
been developed to obtain physical quantities from correlation functions, see |2j S 0> IE S H- The MD approach 
allows the application to large coupling parameters. On the other hand, quantum effects are difficult to include. This 
shortcoming is partially cured by considering pseudopotentials which effectively take into account the uncertainty 
principle by a short distance modification of the Coulomb interaction within the range of the thermal wavelength, see 
More rigorous methods to include quantum effects are wave packet MD simulations (8| or path integral Monte 
Carlo calculations 9J. 

Other points are the finite particle number and the limited accuracy when solving the equations of motion. The 
latter will not be discussed any further. The transition from a finite system to the thermodynamic limit of an infinite 
system can be performed by periodic boundary conditions. The total force on a given particle from all the other 
particles in a basic cell, as well as from the infinite array of their periodic images, can be obtained using the standard 
Ewald procedure 2, ,3] . 

In the present paper, the long-wavelength limit cf{lu) — limfc^o crik^Lu) of the dynamical conductivity is considered 
for a two-component plasma. According to the fluctuation-dissipation theorem (FDT), this transport quantity can 
be expressed in terms of equilibrium correlation functions, in particular the auto-correlation function (ACF) of the 
electrical current or the ACF of the electrical charge density. In the literature see also the internal as 
well as the external conductivity are introduced, relating the electrical current density to the internal or the external 
electrical field strength, respectively. We will present the corresponding relations in the following Section II. An 
important quantity related to the dynamical conductivity is the dynamical collision frequency Analytical 
expressions can be derived in different approximations within a perturbative approach, see |l2j| . 

Section III defines the current ACF in the context of MD simulations, and the connection to the collision frequency 
is shown. While results from MD simulations and analytical approaches for the structure factor and other frequency 
dependent quantities at finite wavenumber k are in good agreement, see e.g. we will discuss the zero-wavenumber 
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case of MD simulations which is relevant for the dielectric function e(fc = 0,uj) or the dynamical conductivity cr(uj). 
Calculations are presented without and including a mean field contribution which lead to the external and internal 
conductivity, respectively. 

Details of the MD simulations are reported in Section IV. Results for the current ACF and the dynamical collision 
frequency at parameter values of a strongly coupled plasma are shown and compared with results of the analytical 
approach. The inclusion of a mean field when performing MD simulations is considered in Section V. The controversy 
between the internal and external conductivity in calculating the collision frequency is resolved. Conclusions are 
drawn in Section VI. 



II. DYNAMICAL CONDUCTIVITY OF THE TWO-COMPONENT PLASMA 



We consider a two-component fully ionized neutral plasma, such as a H plasma consisting of electrons and protons, 
at temperature T and density n of each component. The interaction is given by the Coulomb potential, and the plasma 
is characterised by the nonideality parameter T = e^{4:nne/3y^^{4'!TeokBT)~-^ and the degeneracy parameter Q = 
2mpkj^TTi~'^ {2>i:^n)~^/^ . The linear response to external perturbations in general is presented in various references, 
see e.g. |Tl ll2j . In the following, we will restrict ourselves to relations which are relevant for further discussion. 

Under the influence of an external field Ecxtir,t) = ^gxtgi(fc r-wt) electrical current density (J/c)* is induced. 
The brackets (• • •)* indicate taking the statistical average with the thermodynamic equilibrium distribution. If we 
consider the response in an isotropic systems, the z axis can be selected without loss of generality in such a way that 
Eq'^'^ = Eq^^Cz, k = kcz-, Jk = Jk&z- The relationship between the induced longitudinal current and the external 
field is given by the response function xikj^)- Within linear response theory, x{^t^) is related to the equilibrium 
correlation function of the longitudinal electrical current density 0, 0| 



x(fc,c.) = -^/?^7o-(J^^J^«) (i) 

UJ \ I uj-\-irj 



j die^("+"')*(j^"^(t)J^"s), (2) 



UJ 

where J7o is the normalization volume and the limit rj ^ Q has to be taken after the averaging over the thermodynamic 
equilibrium. Since the longitudinal part of the current density is related to the charge density according to the balance 
equation (due to charge conservation), the longitudinal current ACF can also be expressed in terms of the charge 
density ACF. 

According to the FDT, the response function is related to the dynamical structure factor or the longitudinal part 
of the dielectric tensor e{k,uj) according to (c.f. [iL 13. ITllIT^ 1 

e'°^Hk,u;) ^ 1 - J^^'^'l ^ 1 (3) 
eak^ + X{k,uj) eofc^ 

The longitudinal polarization function Il^°'^^{k,uj) gives the relation between the induced current and the internal 
field as does the dynamical conductivity 

,io„g(fc^^) = ptf°"*^(fc,c.) = . 'f^;. (4) 

It is also called the internal conductivity [lo| . In Eq. H), the dynamical collision frequency v{k,uj) is defined by 
a generalized Drude formula where cijpi ~ (ne^/eoWei)^'^ is the plasma frequency and m^i the reduced mass. The 
phenomenological Drude model is found from the generalized Drude formula Eq. Q if the collision frequency is 
considered to be a real constant equal to the inverse of the relaxation time r in momentum phase space. 

In analogy to the internal conductivity, a so-called external conductivity 'lO] can be introduced from the response 
function ^ 



,{k,cj) = '^x{k,^)^pno{/r^:fr) . (5) 



This quantity is directly related to the longitudinal current ACF. Note that it is not the dynamical conductivity 
defined by Eq. Instead, it is related to the dynamical collision frequency in the following way 

aaxt[k,uj) = — 2— — - . (6) 
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The transverse part of the dielectric tensor can also be related to a conductivity according to 

e(fc,w) = 1 + — a(fc,w). (7) 
The transverse conductivity is defined in analogy to the longitudinal Q as 

a*--(fc,c.) = ^nt--(fc,u;) = . '""j' (8) 

fc^ —tUJ + I/(fc, UJ) 



where i'{k,uj) is commonly called memory function 2]. However, in this case, the Kubo-Greenwood formula |2ll4lllll 
Ek. .14] relates the polarization function directly to the transverse current ACF, 

= pn, ( J*-'^- 4''"")^+,, ■ (9) 

Within a Green function approach, a diagram representation is possible |l2j . In contrast to x{^t^) ^'^d the 
transverse polarization function, which are given by diagrams containing Coulomb interaction in any order, the 
respective current ACF n'"''^^^^^^) 

is given only by the irreducible diagrams. In the long- wavelength limit, transverse 
and longitudinal conductivies lead to the same response of the system, 

lim = lim a'°"«(fc,cj). (10) 



III. CURRENT AUTO-CORRELATION FUNCTION 



Within MD simulations H,^,^^,^, the normalized current ACF 



K{t) = (11) 



is calculated. Here, the long-wavelength limit (fc 0) of the current 

1 ^ 

Jk=,{t) = ^Y.Y.^^<c{t) (12) 

° c i=l 

is considered, where N is the number of electrons and singly ionized ions, and ^ is the speed in z-direction of the 
zth particle of component c, denoted by {i,c}. For convenience, we will drop the index k in the following. Due to 
isotropy, the normalizing factor is equal to 

The Laplace transform of the current ACF reads 

oo 

{J-^J).+^r, = {j') J e^'^-+"^^'K{t)dt. (14) 



On the basis of this quantity, two different results for the conductivity 

oo 

a{uj)^e^u:l,j S''+'^^'K{t)dt (15) 



are derived depending on whether the current densities are considered to be long-wavelength limit of the longitudinal 
or transverse case. 

Firstly, within the transverse response, the Kubo-Greenwood formula © is utilized. The conductivity p5|l is then 
related to the memory function ^{uj) |3, H, HJI via the Drude like formula © and we find 

M= ^o^P' (16) 
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If we assume a constant memory function (collision frequency) v{u!) = the Laplace transformation of a^''^'^^ {u) 
back to K^^'^^^lt) using the functional dependence given by Eq. p5|l . leads to a monotonically decreasing K^'^^^^^t) = 
exp{—vt). This behaviour is observed indeed in simulations for F < 1 0,SIEEj0|- 

Secondly, within longitudinal response, we have to distinguish between the external and the internal conductivity. 
Inserting Eqs. H14|l and into jSJ, this implies that expression (|15|l is the external conductivity. The internal 
conductivity can be calculated via 



^loiip 



O'oxt(t^) 



1 - icroxt(w)/(eow)' 

and due to the generalized Drude formula the collision frequency is, in contrast to H16() . 



LO-pl O-cxt(w) 



Wpi 



UJ 



(17) 



(18) 



Using a constant collision frequency i'{uj) = v in the respective relationship (|SJ| for the external conductivity, we 
find for the longitudinal current ACF via a Laplace transformation 



ifi™^^(t)=exp{-^i} 



2z 



sin(zi) + cos(zt) 



(19) 



This shows that an oscillating behaviour is expected for the ACF. The oscillation frequency tends to Wpi in the limit 
1/ ^ 0. 

If in the long-wavelength limit both v(uj) and v(io) coincide, the current ACF for the longitudinal and transverse 
response cannot be identical. In the following Sections IV,V we will resolve the contradiction between the internal 
conductivity as obtained from the current ACF according to H17|l and (O and the transverse conductivity obtained 
from the current ACF according to H16() and 



IV. SIMULATION TECHNIQUE 



In the MD simulation scheme, the Newtonian equations of motion are solved for a system consisting of iV singly 
charged ions and N electrons exerting Coulomb forces on each other. The ith particle of component c shall be 
denoted as {i,c}. This is a classical treatment where the trajectories of each particle are determined. The original 
Coulomb interaction can be replaced by a pseudopotential, where the short-range part of the interaction is modified 
reflecting the quantum character of the interaction. A systematic derivation of a pseudopotential which reproduces 
the equilibrium properties has been given by Kelbg, see ^ TSj on the basis of the Slater sum. In particular, we use 
the so-called "corrected Kelbg" potential ;15|: 



^cd(r) 



47reo7' 



F 



r 



ScSd 



Acdi^cd) exp - — 




(20) 



where 



1 

rrid ' 



y/2mcdkgT rricd rric 
F{x) = 1 - cxp{-x^) + ^/^x{l - erf(a;)). 



icd 



knTXcd' 



^ee(Cee) = VT^\iee \ + In 



Aeiiiei) = ~\/7rCej + In 



2x/¥|ee 



yexp{-y'^)dy 
exp(7r|^ee|/2;) - 1 



( C(3) + -cmi 



4V^6 



1 - exp(-7r^ei/y) 



where ^(n) are the Riemann-Zeta functions. This interaction potential corresponds to the Coulomb potential at large 
distances and provides the exact value of the Slater sum and its first derivative at r = 0. 

Initially, all the particles are gathered in a cubic box with the edge size L. The number of particles N in this basic 
cell is obtained from a given mean plasma density n via N = nL^. To simulate an infinite homogeneous plasma. 
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images of this charge-neutral basic cell are considered shifting the basic cell by integer multiples of L in different 
directions. This extended system has a constant mean plasma density n. Artefacts may occur due to the periodicity 
of the particle positions, but they are suppressed if the basic cell size is increased. 

The dynamics of both electrons with charge — e, mass mg and ions with charge e, mass rrii is considered. Because 
of the continuous expansion of such plasma, the nearest image method is applied to the force calculation procedure. 



Here, the force Fi 



short 



i long 



on a particle {i, c} is considered to consist of two contributions. The interaction 



forces between particle {i,c} and the nearest neighbour images of all other particles found in the basic cell centered 
around the position r^^c of the considered particle {«, c} is the short-range contribution i^f^™*. The contribution 
is originated from the remaining images, which are not in the basic cell. 
The short-range part of the force is calculated as 



F^ 



short 



N 

EE 



r dV,d{r) 
r dr 



(21) 



The time argument t is suppressed. According to this method it is assumed that the particle {i,c} doesn't interact 
with original particles which at large t may be found far away due to the motion in space beyond the basic cell, 
but with their next neighbours' images obtained by periodically shifting their coordinates into the basic cell centered 
around the particle {i, c}. Thus, the position of each original particle rj^d is replaced by the position of an image 



j,d 



'j,d 



]4 



< 



2' 



(22) 



where a = x,y, z and m is an integer. It should be noted that this procedure is repeated for each particle at fi^c- 
This method implies that each particle is always surrounded by 2A^ — 1 other particles with a constant mean density 
and the plasma is homogeneous in scales larger than the simulation cell. 

The forces F^°'^^ due to the interaction with images outside the basic cell centered around the position ^ of the 
particle {i, c} are treated in a different way. If the dimension L of the basic cell is large in comparison to the screening 
length, the contributions of all images except the nearest one can be neglected. In particular, this is justified in 
the case of a nonideal plasma where the effective interaction potential decreases exponentially with distance due to 
screening. The influence of the far images can be taken into account considering Ewald sums. They are expected 
to give only a small contribution to F^^^°'^^ provided N is high enough. They are not relevant with respect to our 
considerations. 

For explicit MD simulations, we consider a model plasma consisting of singly charged ions and electrons with 
density n = 3.8 x 10^"'^ cm~^ at a temperature of T = 33000 K. This corresponds to recent experiments in dense xenon 
plasmas The plasma parameters introduced in Sec. II take the value T — 1.28, Q — 3.2. It is a nondegenerate, 
strongly coupled plasma. The computations of the current ACF for the ion-electron mass ratios rrii/me = 1836 and 
nii/me = 100 show no considerable difference. Thus the ratio rrii/me = 100 is selected for better convergence when 
averaging over the configurations of ions. The total number of particles N = 250 was found to be enough for F « 1. 
Further increase of the number of particles {N = 400) does not affect any simulation results including the mean 
interaction energy, equilibrium correlation functions and others. The equilibrium state of the plasma at the given 
temperature was obtained using a special procedure described in Q . 

The current ACF is calculated directly from the velocities of the particles in subsequent moments of time according 
to Eqs. Hll|l and H12|l . where fig = with L the length of the basic cell. The averaging of the ACF is performed 
over (1 — 5) • 10^ initial configurations. These configurations are obtained from a long MD trajectory at different time 
moments. As shown in (Q], two configurations are statistically independent if they are taken at times separated by 
the dynamical memory time. In our case about 5 • 10'^ initial configurations are already fully statistically independent 
for electrons. The dynamical memory time for ions increases with the ion mass Thus the smaller mass ratio the 
better averaging for ions is obtained. 

Results are shown in Fig. ^ with circles. The relatively small ion-electron mass ratio {mi/rrie — 100) was chosen 
for computational reasons since the calculation with greater mass ratio shows exactly the same results for the current 
ACF. The current ACF K{t) decreases monotonously as it was also obtained in previous MD simulations 
indicates that the conductivity obtained numerically from K{t) according to H15|) should be treated as the transverse 
conductivity (jS)). The dimensionless dynamical conductivity CT(u;)/(eoa;pi) is shown in Fig. [3 with circles. As w ^ 0, 
the real part has a finite value and the imaginary part vanishes, as expected from Eq. ()16|l . According to the latter 
expression, we then deduct a memory function or collision frequency v(uj) as shown in Fig. |S1 with circles. 

Details of different approximations for the dynamical collision frequency within a generalized linear response theory 
can be found in |12| . The dynamical collision frequency in Born approximation with respect to the statically screened 
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potential (Debye potential) taken in the non-degenerate case and within the long- wavelength limit, is given here 



= 0,w) 



-ign I dy-^ 



V{q) 



WniekBTiloeo 



1 - e^-^'^y 

xy{xy ~ Co — irj) 



(23) 



where 



n 



4^3/2 



Tiu) 



(24) 



In the case of the Fourier transform of the Coulomb interaction V{q) = / [VlQe^q'^) the square brackets become l/y^. 

We will now compare the MD simulations with this analytical treatment of the dynamical collision frequency 
within perturbation theory, see Figs. 0] Firstly, we consider a system with statically screened Coulomb interaction 
V{q) = e^/(rioeo'?^) according to Eq. (|23|) . The results are presented as dotted line. The Born approximation can be 
improved by taking into account the effects of dynamically screening, strong collisions (T matrix) and higher moments 
by introducing a renormalization factor jT3 | in the generalized Drude formula Eq. Q . This approximation is shown 
as solid line. Details of the calculation are given in It can be seen that both real and imaginary part are in good 
agreement with the simulation results for uj < Wpi. This means that in this region the quantum mechanical treatment 
of the Coulomb potential and the classical simulations based on the corrected Kelbg potential are consistent. 

At frequencies w 3> Wpi the asymptotic expansion of the analytical expression for the collision frequency is possible 
using the Fourier transform of the corrected Kelbg potential H20|l 



Vcdiq) 



Seeded 
e^noq 



Erfi(^g)e- 



KdkB 



°^cd(Cc(i) 



qe 



(25) 



For the high frequency behaviour of the real part is found Re i^(w) 



U! 



-3.5 



which is given in Fig. 0] as dashed line. 



There is good agreement between the simulation data and the analytically derived high frequency behaviour. The 
presented analytical treatment was also confirmed by MD calculations of the dynamical structure factor at finite k 
in where the Dcutsch potential was used. 



V. LONGITUDINAL CONDUCTIVITY 



We now investigate the evaluation of the longitudinal conductivity by MD simulations. The current ACF ii''''^"^ 
discussed in the previous Section cannot be taken since this current ACF yields the correct collision frequency only if 
the external conductivity is related to a Drude ansatz. However, this is not consistent. Instead we have to derive the 
internal conductivity, from which a collision frequency can be obtained via the Drude formula. Therefore, the current 
ACF iir'°"s has to be calculated differently than the ACF K^'^'^'^^. It will be shown how to obtain the longitudinal 
current ACF in the long-wavelength limit. However, we note that for finite wavevector k excellent agreement for 
the dynamical structure factor from MD simulation and analytical expressions has been found Q. The condition 
k > 2tt/L means that any charge density wave occurs already within the basic simulation cell and the corresponding 
mean electric field is accurately taken into account. The limit fc is not trivial. For any small k, the system 
is nearly homogeneous, but charge densities (or surface densities) are present at large distances, which can also be 
considered as a mean field. 

For this, we follow the procedure to construct an infinite system by periodic images of a basic cell. We consider 
this as a limiting case of a finite number of images. Denoting the images in z-direction by iVimagcs, then a surface 
of our system is obtained at z- — —{2Nijnases + 1) • L/2 and z+ = (2Afi„iagcs + 1) ■ L/2. When considering the force 
calculation procedure, there are contributions to the forces originating from a surface charge density. This occurs if 
positive and negative charges are moving at different rates across the surface of the basic cell. The introduction of a 
finite number of images compensates this effect at the interfaces, but not at the surface of the whole system including 
all the images. A large dipole moment follows connected with a finite polarization of the system. This surface charge 
density will produce an electrical field which has to be taken into account even in the limit when the number of images 
goes to infinity. If the surface is far away, it produces a homogeneous electrical field E{t) within the simulation box. 
Following this reasoning, it is necessary to include a mean field in the long-wavelength limit as shown below. As a 
consequence, plasma oscillations are obtained in the current ACF. 

On the macroscopic level, the Maxwell equations relate this mean field E{t) to the average current density J{t), 
which is oriented in z-direction according to the conventions in Sec. II, 



(26) 
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Taking the current density according to Eq. H12() as an average over the basic simulation cell and the initial condition 
E{0) = 0, the integration of Eq. leads to 



J3 f-e^^^.e + e^r,,, J 
\ 1=1 1=1 / 



+ (27) 

In this approach, the long-range interaction forces are given by = ecE{t). In particular, the equation of 

motion for an electron includes two parts 

me^ = 4?°"-e^- (28) 

The interaction forces F^^^°'^^ originate from close partners in the Debye sphere within the basic cell. It is fluctu- 
ating around a nearly zero mean value. Nevertheless, the amplitude of these fluctuations are much higher then the 
fluctuations of eE. 

In the MD method, if no mean field term is taken into account, the total energy 

N N N 

= s,., + = ^ E E ^^''(^^..^ - + ^ E <e + f E <^ (29) 

is conserved. If the particle trajectories are calculated including the mean field force, the energy S^^^ + f^jj^ is not 
conserved. Nevertheless, the conservation law can be fulfilled by including the mean field energy £ficid — L^eoE^/2 so 
that the total energy £[^^ — £p^^ + S^.-^^ + f field is conserved. This is illustrated by simulations below. 

The occurrence of plasma oscillations can be demonstrated in the following way. If the mass ratio between electrons 
and ions mi/rrie is large the ion current can be neglected in Eq. p2|l . After that the derivative of the total current 
density is obtained from 

^ = -I5E^ = ;^(^^-^)' (30) 

i—1 

N N N 

^=^E^r*-^EE^.- (31) 

i—1 i—1 j—1 

The force ^ includes only electron-ion interaction forces as all electron-electron interaction forces are compensated 
since they do not change the total momentum of the electrons. Although the force F.^^^™*^ on each electron is typically 
much greater than the force eE from the mean electric field, the average over all electrons is of the same order of 
magnitude as eE. If we now differentiate Eq. p6() and substitute the derivative of the current using H30() , we obtain 
the equation for the mean field 

d^E , ^ 



u^l^E^^l (32) 



On an average, ^ vanishes, so that plasma oscillations are described. The corresponding oscillations in the current 
ACF are obtained from MD simulations as the results below show. 

We now present MD simulations based on the solution of the equations of motion (|28|l in comparison to the MD 
simulations as presented in the previous Sec. IV where the contribution of the mean field —eE was not taken into 
account. The energy conservation is demonstrated in Fig. El according to Eq. H29|l . It can also be seen that the field 
energy ffioid is rather small compared to the particle energy £'^^^. + f^jn- 

Results for the longitudinal and transverse current ACF are shown in Fig.^ After including the mean field into the 
MD simulations, the plasma oscillations in K{t) become well pronounced in contrast to a monotonously decreasing 
behaviour. It should be stressed that the amplitude of these oscillations does not depend on A^. 

The conductivity calculated according to Eq. l(T3|) is shown in Fig. [21 In comparison to the transverse case, the 
conductivity shows a qualitatively different behaviour. The real part following from the MD simulations including 
mean field is zero for zero frequency as is expected from the expression for the external conductivity ©. For the 
case without mean field, Recr has a finite value. In the high frequency limit, both curves coincide. The dynamical 
collision frequencies v{uj) and the memory function ^{lij) calculated from the simulation data for the ACFs are shown 
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in Fig. 13 As pointed out, the results for the Laplace transform of the ACF differ significantly (Figs.[21l- Nevertheless, 
if Eq. (|18|l is used for the collision frequency and Eq. (|16|l for the memory function v{lo) in order to calculate 
the collision frequency, the results for both coincide quite clearly (Fig. |2l . 

Therefore, our analysis showed that the contradiction between the transverse conductivity which should be identical 
with the internal conductivity in the long- wavelength limit and the external conductivity could be resolved if the mean 
field is taken into account. The difference between lu\v{uj) and 1tiiv{uj) in the low frequency limit is caused by the 
numerical error of Ivciviuj) due to substraction of two large terms in Eq. H18|) . 

VI. CONCLUSION 

Molecular dynamics simulations of strongly coupled plasmas were performed using the quasiclassical Kelbg inter- 
action potential. The current auto-correlation function was computed for a non-degenerate two-component plasma. 
Whereas for finite k the dynamical structure factor and the plasma oscillations are reproduced by MD simulations, 
see 0;SI3i the original methods do not allow to consider k values with k < 2ii/L. On the other hand fc = should 
be possible to investigate with MD simulations in a finite volume. 

We presented calculations for the transverse current ACF as well as for the longitudinal one. Although in the limit 
fc — > the transverse and longitudinal dielectric function and conductivities, respectively, coincide, the current ACF 
behave differently in this limiting case. It was shown that the results of MD simulations without a mean field in the 
long wavelength limit provide the monotonously decreasing transverse ACF. Its Laplace transform is to be directly 
related to the transversal conductivity. 

In MD simulations for the longitudinal case, a mean-field term has to be included into the equations of motion in 
addition to the short range forces inside the Debye sphere. This mean-field term originates from surface charges not 
taken into account in the usual procedure of force calculation by the nearest image method. Simulations with these 
altered equations of motion show well pronounced plasma oscillations in the longitudinal current ACF. The results 
for the collision frequency as obtained in both simulation methods using the corresponding relations for the internal 
or external conductivities do coincide. 

Additionally, the dynamical collision frequency inferred from the simulation data was compared with analytical 
results, which were derived using a generalized linear response theory. We found good agreement in the low and high 
frequency limits for a moderate nonideality. In particular, for w < Wpi, classical MD simulations using the corrected 
Kelbg potential are able to reproduce the quantum behaviour of Coulomb plasmas. 
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FIG. 2: Real and imaginary parts of tlic Laplace transformation of the current ACF; MD simulations without (circles) and 
including (triangles) an additional mean-field term in the equations of motion. 
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FIG. 3: Real and imaginary parts of the dynamic collision frequency or memory function from MD simulations without (circles) 
and including (triangles) an additional mean-field term in the equations of motion. 
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FIG. 4: Dynamical collision frequency within different methods; points - MD simulations; analytical approximations: dotted 
line - Born approximation Eq. II23II with Coulomb potential, solid line - same approach including dynamically screening and 
strong collisions (T matrix) and higher moments via renormalization factor IT3|, dashed line - high frequency asymptote for 
Born approximation Eq. 12311 with corrected Kelbg potential. 
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FIG. 5: Conservation of the total energy in MD simulations; curve 1 - total energy of the particles Ep^t + according to 
Eq. 12911 . curve 2 - total energy El^t including the mean field energy. 



